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Motivated by the discovery of gravitational waves (GWs) at aLIGO and no evidence of new physics 
at current LHC, we discuss that a generic classes of extended new physics models with hidden gauge 
group could undergo one or several times first-order phase transitions associated with the gauge 
group symmetry breaking during the evolution of the universe, which might produce detectable 
phase transition GWs signals at future GWs experiments, such as eLISA and BBO. 
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I. INTRODUCTION 


The observation of gravitational waves (GWs) by Advanced Laser Interferometer Gravitational Wave Observatory 
(aLIGO) has initiated a new era of exploring the cosmology, the nature of gravity as well as the fundamental 
particle physics by the GWs detectors [2H8]. Especially, due to the limitation of the colliders’ energy, GWs detectors 
can be used as new techniques to probe the symmetry breaking patterns or phase transition history for a large class 
of new physics models with hidden gauge group extensions of the standard model (SM), which are motivated by 
the mysterious experimental results in our understanding of particle cosmology (such as the dark matter problem or 
the puzzling observed baryon asymmetry of the universe), and the absence of new physics signals at current collider 
experiments. The increasingly attractive new physics models with hidden gauge group extensions of the SM have many 
new particles without leaving any observable imprints at current particle colliders. However, the GWs experiments may 
provide a possible approach to test their existence. For example, to explain the baryon asymmetry of the universe via 
electroweak baryogenesis by the hidden gauge group extended model, a strong first-order phase transition is needed to 
realize the departure from thermal equilibrium [9HIT]. And during the first-order phase transition, the detectable GWs 
will be produced through three mechanisms: collisions of expanding bubbles, sounds waves, and magnetohydrodynamic 
turbulence of bubbles in the hot plasma 19|. Phase transitions in particle physics and cosmology are usually 
associated with the symmetry breaking, i.e. where the universe transits from a symmetric phase to a symmetry broken 
phase when the temperature drops below the corresponding critical temperature. 

For the first time, we have a chance to explore the hidden gauge symmetry breaking through phase transition GWs 
signals after the discovery of the GWs by aLIGO, which is particularly exciting. In this paper, we study the possibility 
to explore the hidden gauge symmetry breaking patterns by the phase transition GWs signals. In particular, we will 
focalize our analysis to GWs detection of the non-Abelian gauge group extended models, where the symmetry breaking 
at each energy scale may associate with a first-order phase transition, as shown in Fig[I] The hidden group Gyidden; 
which can spontaneously breaks into the SM gauge group Gsm : SU(3), ® SU(2); ® U(1)x. During the hidden gauge 
symmetry breaking, a strong first-order phase transition can take place, which will produce detectable phase transition 
GWs. For example, the hidden gauge group G'Hidden can be the non-Abelian group SU(3), ® SU(3); ® U(1)y, and 
the model is called 3-3-1 model [21]. We show that many versions of the 3-3-1 model can produce at least one 
time strong first-order phase transitions at the energy of several TeV in some parameter spaces, which can produce 
detectable GWs spectrum by Evolved Laser Interferometer Space Antenna (eLISA) [22], Big Bang Observer (BBO) [23]. 
In general, there can exist several times spontaneous symmetry breaking, which may also accompany several times 
first-order phase transition during the evolution of the universe as shown in Fig]1} If the energy scale of the symmetry 
breaking and the associated first-order phase transition range from 107 to 10° GeV triggered by the hidden sector, 
the phase transition GWs spectrum may be within the sensitivity of future aLIGO. Other attractive examples of 
GWs detection of the non-Abelian gauge group extended models are also discussed, such as the first-order dark QCD 
phase transition GWs in the relaxion mechanism [24], which may be tested by pulsa time array (PTA) at the Square 
Kilometre Array (SKA) [25] or the Five-hundred-meter Aperture Spherical Telescope (FAST) [26] in China. The other 
one is the Nnaturalness mechanism [27], which may also produce observable phase transition GWs signals at the GWs 
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FIG. 1: Symmetry breaking (phase transition) patterns for the non-Abelian gauge group extended models with the evolution of 
our universe, where the first-order phase transitions may happen. 


detector. 

This paper is organized as follows: In Section II, we schematically discuss the GWs detection of the hidden gauge 
group symmetry breaking and show how to calculate the phase transition GWs during the first-order phase transition. 
In Section III, we will study the phase transition GWs spectra in some concrete hidden gauge group extended models. 
In Section IV, we show our final discussions and conclusions. 


II. FIRST-ORDER PHASE TRANSITION GWS SPECTRUM 


In a large class of new physics models with extended non-Abelian hidden gauge group, one or several times 
cosmological phase transitions can occur during each step’s hidden gauge symmetry breaking. In general, we can have 
the gauge symmetry breaking patterns with the extended hidden gauge group at each scales with the evolution of the 
universe as shown in Fig.|1| For example, the hidden symmetry breaking pattern may be the following: 


G(Hidden N) -.- > G(Hidden1) > G(SU(3)c@SU(3)L@U(1)y) > G(SU(3)c@SU(2)L@U(1)x) > G(SU(3)c@U(1)zm). 


With the evolution of our universe, symmetry breaking will happen at corresponding energy scale, where the strong 
first-order phase transition may take place. Detailed models are given in Section III. 

Along with the symmetry breaking from a ‘false’ vacuum to a ‘true’ vacuum, strong first-order phase transition 
occurs if there exists a sufficient potential barrier between them. This process can produce an observable stochastic 
GWs signals, which can be detected in some GWs detectors, such as aLIGO, eLISA, BBO, SKA, FAST and so on. 
Their sensitivity range for some critical temperatures depends on the energy scale of the first-order phase transition 
for different hidden gauge group extended models, as shown in section III. 

To discuss the GWs spectrum from the first-order phase transition, it is necessary to begin with the one-loop finite 
temperature effective potential: 


Ver (®, T) = Viree(®) + Vew(®) + Vihermal(®, T) + Vaaisy (®, T), (1) 


where ® represents the order parameter of the phase transition (a real scalar field), Vow(®) is the one-loop Coleman- 
Weinberg potential at T = 0, and Vihermai(®,T) + Vaaisy(®, T) is the thermal contribution including the daisy 
resummation [28]. 

During each step of symmetry breaking in the hidden gauge group extended models, strong first-order phase 
transition may occur. During the first-order phase transition, bubbles are nucleated via quantum tunneling or thermally 
fluctuating the potential barrier. The bubble nucleation rate per unit volume T is given by I =Io(T)e~ 9") with 
To(T) x T4 [29], where Sz(T) is the Euclidean action [80] BI] defined as 


se(t) = | drëx E (=) + 5(V®)? + Vor(®,7) (2) 


Here, Sp(T) ~ $3(T)/T, and T = Tpe—3/T [29] where 


s(t) = f ax | 5(V9)? + VaT). (3) 


From the above equations, in order to obtain the bubble nucleation rate, the profile of the scalar field ® during the 
bubble nucleation needs to be calculated by solving the following bounce equation: 


Pb 2d VælÈ,T) 


T = 4 
dr? rdr oğ i (4) 
with the boundary conditions 

d® 

S (r=0)=0, (5) 

(r = œ) = Praise. (6) 


The bounce equation can be solved numerically using the overshoot/undershoot method. The first-order phase 
transition terminates when nucleation probability of one bubble per horizon volume is of O(1), i.e., '(T,) ~ H4. That 
is to say, it should satisfy 


53(T.)/T. = 41n(T, /100GeV) + 137. (7) 


It is known that there exist three major sources for producing GWs during the first-order phase transition, which 
are collisions of the vacuum bubbles [I5], sound waves and turbulence [I7] [78] in the plasma after collisions, 
respectively. 

The most well-known source is the bubbles collisions, and the first-order phase transition GWs spectrum depends on 
four parameters. The first one is the ratio a of the vacuum energy density released in the phase transition to that of 
the thermal bath, defined as 


T; 

oa =, (8) 
Prad (Ta) 

where * specifies that the quantity is evaluated at the time Tą determined by Eq. (7). Note that, the false vacuum 


dv? 
2 aT . 
Praa (Tx) is equal to 359+(T)T*. Here, eff (T) is the temperature-dependent true minimum of the effective potential 
of the scalar field, which is responsible for the phase transition. The parameter a measures the strength of the phase 
transition GWs, namely, larger values for œ correspond to stronger first-order phase transition GWs. 
The second one is the time duration of the phase transition 871, where one has 8 = — Be | on is \:<z,- In other 


energy density ¢(T,) = [T —Vmn(T)||r=r, is the latent heat density, and the plasma thermal energy density 


words, 871 corresponds to the typical time scale of the phase transition. Since 6 = È/ T during the phase transition 
from its definition, one has 


B _ a($3/T) 
«(aE 


(9) 


T=T, 


The third one is the efficiency factor A, which characterizes the fraction of the energy density converted into the 
motion of the colliding bubble walls and the last one is the bubble wall velocity vy. The energy released into the GWs 
of peak frequency is then given by [32]: 


2 2 
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The second and third sources are the GWs from the matter fluid effects, which can further contribute to the total 
energy released in gravitational radiation during the phase transition. Here, we just use the formulae given in [33]. One 
source is from the sound waves in the fluid, where a certain fraction A, of the bubble wall energy (after the collision) is 
converted into motion of the fluid (and is only later dissipated) with the following contribution 


PGW,sw H, 2 a? 
Z w Osw | — ] A, | ——; | - 11 
Ptot : (F) Ta m 


The other one is from turbulence in the fluid, where a certain fraction Atu of the walls energy is converted into 


turbulence, with the contribution 
PGW,turb H. 3/2 o3/2 
< w Oty, | — ) A —,, }. 12 
Ptot : (5) ™ (a a 


It is worth noticing that these two contributions from the matter fluid effects depend on H,,./ linearly, and they are 
not fully understood. In some cases, these two effects may be larger than the one from bubble collisions. 

The peak frequency produced by bubble wall collisions at T, during the first-order order phase transition is given 
by [36] [37]: fz, = 0.628/(1.8 — 0.1u, + v2). Considering the adiabatic expansion of our universe from the early universe 
to the present universe, the ratio of scale factors at the time of first-order phase transition and today can be written as 


a sO K0 Haa i (sacar) (=) (13) 


where gt is the effective degrees of freedom at Tą. Thus, the peak frequency today is feo = f*,ax/ao, and the 
corresponding GWs intensity is given by 


Oo (fh? ~1.67 x 10-* (2) "( da Wey 
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The peak frequency of the GWs signals due to the sound wave effects is about f*, = 28/(V3v») [16 B3], and 
similarly its current value is fsw = fx,@x/ao. In this case, the GWs intensity is expressed as [33] 


yy (f)h? 2.65 x 10-6 (FE (Ave)? (100), 
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in which A, ~ a (0.73 + 0.083/a + a) ' for relativistic bubbles [34]. 

The GWs signals produced from the turbulence in the plasma have the peak frequency at about fă, = 1.758/vp 
[33], which determines the current value as ftu = f%,a«/ao after considering the red shift. The phase transition GWs 
intensity from turbulence is formulated by 


Qa ( fh? =3.35 x 10-4(=) (wey) oi 


` (F/F)? 
OF Ffa) O + 8a Fao/(a.H)) 


The final phase transition spectrum consists of the three contributions above. 


II. PHASE TRANSITION GWS FROM NON-ABELIAN GAUGE GROUP EXTENDED MODELS 


In this section, we discuss the phase transition GWs in some non-Abelian gauge group extended models in details, 
where one or several times symmetry breaking and first-order phase transition may occur with the evolution of our 
universe at certain critical temperature. Firstly, the GWs spectra in the hidden gauge group extended models based 
on the SU(3), ® SU(3), Q U(1)y gauge symmetry, commonly known as the 3-3-1 model [20] [21] are investigated. 
The 3-3-1 models can explain the electric charge quantization and three generations of fermions [20] BI]. The collider 
phenomenology of the 3-3-1 models has been extensively studied, such as the recent Ref. and references therein, 
and the phase transitions in some versions of 3-3-1 models have been studied by Refs. [39141]. Here, we try to explore 
the possible phase transition patterns associated with the symmetry breaking SU(3), ® U(1)y > SU(2); @U(1)x > 
U(1)em using the phase transition GWs in some versions of the 3-3-1 models, where the scalars fields are accommodated 
in a certain representation of the SU(3); gauge group in each version. 

In the following, we will show three versions of 3-3-1 models (the minimal, the economical and the reduced minimal 
3-3-1 models, respectively) [39}{41], where very strong first-order phase transitions can take place at the TeV scale [39} 
AI] and produce detectable phase transition GWs by eLISA and BBO. For simplicity, we limit our discussion of the 
first-order phase transition to the thermal barrier case, where the potential barrier in the finite temperature effective 
potential origins from thermal loop effects. In this type of thermal barrier case, the bosonic fields contribute to the 


thermal effective potential of the form Veg Ə (—T'/12m) (me son(X, D * in the limit of high-temperature expansion. 


To begin with the concrete prediction of the GWs in the following examples, we first discuss the general effective 
potential near the symmetry breaking and phase transition temperature, which can be approximated by 


1 eT À 
VælX,T) ~ 5 (—p? + eT?) X? - ee a + re (14) 
Here, X represents the field of the order parameter for the phase transition. For the electroweak phase transition 
in the SM, X field is just the Higgs field. The parameter e quantify the interactions between X field and the light 
bosons, and can be schematically written as e ~ ) right bosonic fields (degrees of freedom) x (coupling to X)*/?. And, the 
parameters c depends on interaction between X and light bosons and fermions. For the heavy fields whose masses are 
much larger than the critical temperature, their contribution can be omitted from Boltzmann suppression. This can 
help to simplify our discussions when the models have many new fields at different energy scales. Thus, in this case of 
qualitative analysis, the wash out parameter can be obtained as 


(X(T) e 


T. x 67. ? (15) 


where the angle bracket means the vacuum expectation value (VEV) of the field X at the critical temperature Te. 
From the above qualitative analysis, we know that introducing new light bosonic fields (compared to the corresponding 
critical temperature) helps to produce or enhance the first-order phase transition. The 3-3-1 models just introduce 
enough bosonic fields to produce detectable phase transition GWs. 


A. Phase transition GWs spectrum in the economical 3-3-1 model and the reduced minimal 3-3-1 model 


We consider the first-order phase transition GWs spectrum in the so-called economical 3-3-1 model [40]. In this 
version of 3-3-1 model, one chooses the simplest SU (3)z representations for the scalar fields with spontaneously 
symmetry breaking, namely, two complex scalar triplets with different hypercharge are needed as follows 


1 
X= (x9,x9 x8)" ~ (3-3) , (16) 


T 2 
b= (0t) ~ (33). (17) 
The scalar potential is written as 


V(X% e) = ÉX XHA? + ug e + Algo) 
+A OOl e) + Aaltola. (18) 


The SU (3) ® U(1)y gauge group is broken spontaneously via two steps. In the first step, the symmetry breaking 
SU (3) 9 U(1)y > SU (2) ® U(1)x happens when the scalar triplet x acquired with the VEV given by 


1 T 
= — (u, 0,w) . 19 
(xX) i ( ) (19) 
In the last step, to break into the SM U(1)em gauge group SU(2); S U(1)x — U(1)gm, another scalar triplet ¢ is 
needed to acquire the VEV as follows 


(8) = 0.0.0)", (20) 


In this version of 3-3-1 model, there exist two neutral scalars!, one is the SM Higgs boson h, the other is the heavy 
scalar Hı. In this paper, the packet ‘CosmoTransitions’ [42] is used to numerically calculate the first-order phase 
transition. During the first symmetry breaking SU (3)z @ U(1)y > SU(2)z ® U(1)x, the order parameter field for the 


1 Here, we do not present the details of the other interactions and the whole particle spectra. Details on the parameter space discussions of 
phase transition in this model can see Ref. [40]. 
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FIG. 2: The phase transition GWs spectra h?Qaw for the economical 3-3-1 model. The colored regions correspond to the 
expected sensitivities of GWs interferometers eLISA and BBO, respectively. The black line depicts the GWs spectra for 
the two times first-order phase transition during SU (3)z @ U(1)y = > SU(2)z ® U(1)x at the TeV (right line) scale and 
SU(2)z Q U(1)x = > U(1)gm at the electroweak scale (left line), respectively. 


phase transition is approximately the H; scalar field, namely, the X = H; field if it is compared to Eq. (14). Strong 
first-order phase transition at the energy scale of several TeV can be induced by the new bosons and exotic quarks in 
this version of 3-3-1 model if the mass of these new particle are from 10? GeV to several TeV. During the last symmetry 
breaking SU(2), ® U(1)x — U(1)gm, the order parameter field for the phase transition is approximately the Higgs 
field, namely, the X = h field if it is compared to Eq. (14). First-order phase transition at the electroweak scale can be 
triggered by the new bosons and it is rather weak than the first step. For the set of benchmark points? of my, = 1.4 
TeV, m H= 3.2 TeV, using the methods and formulae in the above Section |II|and the package ‘CosmoTransitions’ , 
these two times first-order phase transition will produce two copies GWs spectra with different characteristic peak 
frequency and amplitude, as shown in Fig.|2} It shows that the phase transition GWs can be used to explore the 
symmetry breaking and phase transition patterns by eLISA and BBO. 

The GWs spectra of the reduced minimal 3-3-1 model is similar to the one of the economical 3-3-1 model since their 
symmetry breaking and phase transition patterns are similar to the economical model. The reduced minimal 3-3-1 
model is mainly composed by neutral scalars h, Hı, doubly charged scalar ht*, two SM like bosons Z1, WF, the new 
heavy neutral boson Z2, the singly and doubly charged boson U** and USS. These new particles and exotic quarks 
can be triggers for the first-order phase transition [41]. Taking use of the methods and formulae in See[I]] and the 
package ‘CosmoTransitions’ [42], for the set of benchmark points my, = 1.3 TeV, mp++ = 3.3 TeV, two copies GWs 
spectra during two times first-order phase transition will be produced, as shown in Fig.|3} which can also be detected 
by eLISA and BBO. 


B. Phase transition GWs spectrum in the minimal 3-3-1 model 


The minimal 3-3-1 model [39] corresponds to the parameter 6 = —\/3 [89] in the definition of the electric charge 
operator. The gauge bosons, associated with the gauge symmetry SU(3), ® U(1)y of the model, consist of an octet 
wi (i =1,--- ,8) and a singlet B„. In this version of 3-3-1 model, three SU (3)z triplets scalars are needed to break 
the symmetry and generate the mass of gauge bosons and exotic quarks as the following 


= (x7 x x)". (21) 


2 We show only one set of benchmark points to avoid confusion between the GWs spectra from two times first-order phase transition and 
the GWs spectra for one time first-order phase transition with two sets of benchmark points. 
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FIG. 3: The GWs spectra h?Qaw for the reduced minimal 3-3-1 model. The colored regions correspond to the expected 
sensitivities of GW interferometers eLISA and BBO, respectively. The black line depicts the GWs spectrum for the two times first- 
order phase transition during SU(3), QU (1)y = > SU(2)L 8U (1)x at the TeV scale (right line) and SU (2) ®U(1)x = > U(1)em 
at the electroweak scale (left line), respectively. 


The scalar potential for p, 7 and x is [39] [44] 


V (p,m, x) = panty +A (nin)? + p3ptp 
+ rz (ptp)” + paxty +As (xtx)” 


+ o (øp) (xx) + Az (o 
+ As xin) (ntx) +o (etx) htp) 


1 B 
Po (fie nipjxk +H. c.) ‘ (22) 


The new gauge bosons acquire masses at several TeV scale when the SU(3)z x U(1)y group breaks down to 
SU(2)z x U(1)x triggered by the SU(3)z scalar triplet x, while the ordinary quarks and SM gauge bosons obtain 
their masses at the last step symmetry breaking trigged by the triplets 7 and p. There exist three neutral scalars 
including the lightest one which corresponds to the SM Higgs boson h and the other heavier scalar bosons H? and H9. 
By calculating the one-loop effective potential based on Eq. (1) and taking use of the package ‘CosmoTransitions’ [42] 
in this minimal 3-3-1 model, we find that there are regions of parameters allowed by the collider constraints that 
can favor a strong first-order phase transition, when the gauge group spontaneously breaks from SU(3), x U(1)y to 
SU(2)z x U(1)x as shown in Ref. [39]. During this phase transition, the X field here is just the H? field. Then, the 
phase transition GWs spectrum can be obtained from the above GWs spectrum formulae. Since this model has so 
many free parameters, which makes it very complicated to study the whole parameter regions allowed, we only show 
one set of benchmark points. For example, when vy, = 4 TeV and Mpo = 0.8 TeV, the corresponding GWs spectrum 
is shown in Fig. |4| which is within the sensitivities of eLISA and BBO. 


C. Phase transition GWs spectrum in other hidden gauge group extended models 


In general, if the SM is extended by hidden non-Abelian gauge group, strong first-order phase transition may 
occur associated with the symmetry breaking at each symmetry breaking scale, where the phase transition GWs may 
be produced. Then, phase transition GWs can be used to test the hidden gauge symmetry breaking. Besides the 
above phase transition GWs in the 3-3-1 models, we first discuss the new mechanism, namely, the so called relaxion 
mechanism [24] proposed in 2015, where the light Higgs masss comes from the dynamical cosmological evolution during 
the early universe. The relaxion mechanism can technically relax the EW hierarch and the highest cutoff 
relaxed in Ref. is about 10° GeV. The original relaxion mechanism in Ref. need the inflation sector and the 
QCD.-like sector. Especially, one solution to avoid the strong CP problem in the simplest relaxion model, the dark 
QCD gauge group is needed, where new dark fermions are also included. In some allowed parameter space, the dark 
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FIG. 4: The GWs spectra h? Qaw for the minimal 3-3-1 model. The colored regions correspond to the expected sensitivities of 
GWs interferometers eLISA and BBO, respectively. The black line depicts the GWs spectrum for the first-order phase transition 
during SU (3)z ® U(1)y = > SU(2)r ® U(1)x at the TeV scale. 


SU(3)dark Can give first-order QCD phase transition at the dark QCD scale. If the dark QCD scale is about O(100) 
MeV, the first-order phase transition can produce phase transition GWs with the peak frequency in the 107° — 1077 
Hz range, which can be probed by the PTA GWs experiments, such as the SKA or FAST®. A schematic GWs spectrum 
for the dark QCD phase transition is shown in Fig.|5| Very recently, a novel mechanism called “Nnaturalness” [27] 
has been proposed to solve the electroweak hierarchy problem by introducing N copies of the SM with varying values 
of the Higgs boson mass parameter in a unique universe. Although, the special reheaton particle has been added to 
the “N” naturalness mechanism to suppress the baryogenesis from other copies with v 4 246 GeV, the first-order phase 
transition may still happen in some parameter spaces and produce detectable GWs signals, which will be carefully 
discussed in our future work [46]. However, we show that if a first-order phase transition takes place at a critical 
temperature of O(107—10°) GeV [A [47], such as some versions of grand unified models, this could potentially produce 
detectable GWs spectrum in the future aLIGO experiments, and provide us with a unique probe of the hidden gauge 
symmetry breaking at high energy scales. 


3 The phase transition GWs signals in dark matter models with SU(N) hidden gauge group are discussed in By. 
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FIG. 5: The phase transition GWs spectra h?Qaw from first-order phase transition during the evolution of our universe. The 
colored regions represent the expected sensitivities of GWs detectors aLIGO, eLISA, BBO and SKA, respectively. The red line 
depicts the possible GWs spectrum where the first-order phase transition occurs at the energy scale of O(100) MeV in some 
dark QCD models. The black line represents the GWs spectrum for the first-order phase transition at about 4 TeV in some 
dark gauge group models, such as some 3-3-1 models. The purple line corresponds to the GWs spectrum when the first-order 
phase transition takes place at the energy scale of O(10000) TeV in some models with high symmetry breaking scale. 


IV. DISCUSSIONS AND CONCLUSIONS 


In Fig.|5} the first-order phase transition GWs spectra h?Qqw during the evolution of our universe are shown for a 
generic classes of hidden gauge group extended models at the different energy scales. The colored regions represent the 
expected sensitivities of GWs detectors SKA, BBO, eLISA and aLIGO, respectively. The red line depicts the possible 
GWs spectrum where the first-order phase transition occurs at the energy scale of O(100) MeV in some dark QCD 
models. The GWs produced in the dark QCD models can be detected by the PTA experiments, such as the planned 
SKA or the FAST built in China. The black line represents the GWs spectrum for the first-order phase transition at 
about 4 TeV in some dark gauge group extended models, such as some versions of 3-3-1 models. We have shown that 
all the three versions of the 3-3-1 models discussed above could produce strong first-order phase transition GWs at 
TeV scale when the dark gauge symmetry SU (3)z ® U(1)y breaks to SU(2);, @ U(1)x. Especially, in the economical 
and reduced minimal 3-3-1 models two times first-order phase transition may take place, which will produce two GWs 
spectra with different characteristic peak frequencies. In general, the phase transition GWs produced at the scale 
from O(100) GeV to several TeV can be tested at future laser interferometer GWs detectors in space, such as the 
planned eLISA, BBO, Taiji and Tianqin [48]. The purple line corresponds to the GWs spectrum when the first-order 
phase transition takes place at the energy scale of O(10000) TeV in some high scale models, which may be within 
the sensitivity of the future aLIGO and provide us with a unique detection of the hidden gauge symmetry breaking 
at high energy scales beyond the abilities of particle colliders. Besides these schematic models, many other hidden 
gauge group extended models may also undergo one or several times first-order phase transition at different energy 
scale as shown in Fig[i] where the corresponding phase transition GWs spectra can be produced and tested at the 
corresponding GWs detectors. 

To conclude, phase transition GWs becomes a new and realistic approach to explore the the particle cosmology and 
fundamental physics. From the aspect of cosmology, our universe may undergo one or several times phase transition 
during the early evolution of the universe. And we can hear the cosmological phase transition using the phase transition 
GWs if there exists first-order phase transition. From the aspect of particle physics, this phase transition GWs 
approach can compensate for the collider experiments to explore the hidden gauge group extended models and provide 
a novel approach to probe the symmetry breaking or phase transition patterns. More detailed study will be discussed 
in the future. 
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